Probing for the Charm Content of B and T Mesons 
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A slow J/ip bump exists in the inclusive B — > J/ip + X spectrum, while the softness of J/ip spectrum 
in Y(IS') — » J/ip + X decay is in strong contrast with expectations from color octet mechanism. We 
propose intrinsic charm as the explanation: the former is due to B — > J/ipDir, with three charm 
quarks in the final state; the latter is just a small fraction of T(15) — > (cc) s i ow + 2 "jet" events, 
where the slow moving cc system evolves into pairs. Experimental search for these phenomena 
at B Factories and the Tevatron is strongly urged, as the implications go beyond QCD. 
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Owing to its heaviness and narrow width, the J / ip me- 
son has helped shed much light on the underpinnings of 
Quantum Chromodynamics (QCD), the more recent ex- 
ample being the inclusive production of energetic J/ip's 
in various processes. Thus, the inclusive spectrum of 
B —>■ J/ip + X decay is largely understood, except for 
a slow J/ip bump observed by CLEO and recently 
confirmed Q by the Belle Collaboration (see Fig. 1(a)). 
Another intriguing old result from CLEO could be re- 
lated, namely the peaking of pj/^ below 2 GeV for in- 
clusive T(IS') — > J/ip + X decay j|, which is in strong 
contrast (see Fig. 1(b)) with the hard spectrum expected 
from color octet mechanism. In this Letter we attempt 
at linking these two soft J/ip production phenomena and 
propose specific, testable mechanisms. 

We suggest that intrinsic charm (IC) content may be 
better probed in heavy hadrons such as B mesons. Intrin- 
sic, in contrast to extrinsic, charm content of hadrons Q 
carry large momentum fraction and are in the lowest 
energy configurations as they arise from energy fluctu- 
ations. Thus, in the decay of heavy mesons, the produc- 
tion of slow J/ip and r\ c mesons is favored. From the indi- 
cation that IC in proton could be ~ 1% ||, we find that 
IC in B could lead to B — > J/ipDn (i.e. ccc in B decay 
final state f) at the few x 10~ 4 level and can in principle 
account for the above excess. FromT(lS') — > J/ip\ so f t +X 
we argue that T(1S) -> £>W| slow + X could be quite 
abundant. These modes are rather promising at the B 
Factories (the T(1S) especially at CLEO) and at the 
Tevatron, the observation of which would provide strong 
impetus for establishing intrinsic charm in hadrons. In 
turn, this would have implications for the flavor program, 
such as the extraction of V c b from B —> D*£v. 

Using 1.12 fb _1 data on T(4S) resonance, CLEO has 
published the inclusive J/ip spectrum from B decays . 
The result, with feed-down from B — > ip' + X (ip' — > 
J/ipirir) and \ c + X (x c ~ > J/ipj) subtracted, is repro- 
duced as stars with errors in Fig. 1(a). Evidently there is 
some activity around Pj/^ ~ 0.4-0.6 GeV. In 2000, Belle 
gave preliminary results H on inclusive B — > J/ip + X 



based on 6.2 fb _1 data. Subtracting the B — > ip' +X and 
Xc + X feed-down by ourselves, we plot in Fig. 1(a) the 
result as diamonds without errors, which confirms that 
there is some activity below 0.8 GeV. 

To bring home the point, some modeling of the inclu- 
sive spectrum has to be taken. For simplicity, we adopt 
the "modified phase space" approach |6) and modulate a 
constant matrix element squared by 



f(p)=p(p^-p)/pl^ e'^ 2 /^ 



(1) 



where p n 



1.95 GeV is the maximum J/ip momen- 
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FIG. 1. (a) Inclusive B — > J/ip + X spectrum with 
B —> ip' +X and Xc + X feed-down subtracted. The stars (di- 
amonds) with(out) errors are from Ref. [1] ([2]). The dashed 
curve is a simple modified phase space fit. The solid line fits 
below D* threshold assuming recoiling Dir; D recoil is also 
indicated, (b) Inclusive T(IS') — > J/ip + X spectrum from 
CLEO [3]. The curve is from color octet mechanism. 
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turn, and po and <7o are adjustable parameters. Taking 
a simple average (not plotted) of CLEO and Belle data, 
we adjust p , <r to 1.9, 0.8 GeV and give the dashed line 
in Fig. 1(a) as a plausible fit. The apparent excess in 0.3 
GeV < Pj/^ i$ 0.8 GeV is of order 5 x 10~ 4 , comparable 
to the rate for B — » J/ifjKs- 

We stress that more sophisticated models do not 
change the above result. For example, a Fermi motion 
model with pp ~ 0.57 GeV for b quark inside the B 
meson can give a good fit to data above 0.8 GeV. 
A parton-based model with soft b quark momentum can 
also work, albeit less well. In both cases one invokes 
NRQCD and the color octet mechanism, but one is un- 
able to fit the low pj/$ excess. In fact, these models have 
a softer lowpj/^ tail than the simple approach of Eq. (1), 
making the excess even more striking. 

The D and D* recoil thresholds at Pj/$ = 0.88 and 
0.66 GeV are indicated in Fig. 1(a). There is no ap- 
parent excess for B — > J/ipD. Taking into account the 
broadening from B motion in T(4S*) frame, Belle data 
might indicate the presence of B — * J/tpD* ~ 10~ 4 . 
Since D* marks the opening of the Dir threshold, we 
adapt Eq. (1) to p max , Po, Co = 0.66, 1.4, 1.0 GeV, and 
fit with the solid curve in Fig. 1(a). It appears that 
B — > J/ipDir ~4x 10~ 4 could account for the lump at 
low pj/if, With three charm quarks in the final state, 
it would be rather distinct and should be searched for. 
The challenge, however, is to account for such rates. 

Two possible diagrams are given in Figs. 2(a) and 
2(b). The first involves W exchange and is of annihi- 
lation type. Since nonperturbative cc production is ex- 
ponentially suppressed, the leading mechanism is pertur- 
bative and via one gluon. Collecting factors, we esti- 
mate that the rate should be suppressed by (/s/ms) 2 x 
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FIG. 2. Standard B° -> J/?/;D ( * )0 via (a) exchange 
and (b) rescattering; B — > J/ipD^*' (jfj via (c) exchange and 
(d) spectator from \bccq) Fock component, and (e) be annihi- 
lation; (f) T(1S) — * J/ip + X via \bccb) Fock component. 



(a s /n) 2 x 1/3 compared to B — > DD S ~ 1%, where 
Pj/ifj/pD 3 ~ 1/3 comes from phase space. This gives 
B" -> J/ipDM° ~ 10~ 7 from Fig. 2(a) alone. 

Fig. 2(b) is much harder to estimate, since it involves 
dd — ► cc rescattering. Assuming that this occurs per- 
turbatively via one gluon, we estimate the rate from 
Fig. 2(b) alone to be no more than (j x (a s /Tr) 2 x 1/3 
times B — » DD~ rate ~ 1%. Since the rescattering de- 
mands energetic d quark and spectator d quark to pro- 
duce the heavy cc system, the factor is expected to 
be considerably less than one, hence B° — > J/ipD^ 
from Fig. 2(b) should be safely below 10 -5 . As this 
is not a firm result, we note that if one replaces cc by 
ss in Figs. 2(a) and 2(b), one already has the limit 
B° -> D+K- < 10~ 4 from data g. This can be viewed 
as an upper bound on B° — > J/ipD^ modes, where 
the limit can be improved at the B Factories. Similarly, 
rescatterings such as B° — > D + ir~ — > J/tpD (corre- 
sponding to two cuts of Fig. 2(b)) should not only suffer 
from cancellations between the large number of possible 
diagrams, but can also be firmly bound by searching for 
color suppressed modes such as D°ir°. 

We thus see that, while hard predictions of hadronic 
B decays are difficult, there would likely be no plausible 
explanation for B — » J/i/jDtt if it is observed at xlO -4 
level or higher. 

It is intriguing that a soft J/tp problem exists for T(15*) 
decays as well. Based on ~ 7x 10 5 T(15) events and ~ 20 
J/tp — > candidates, CLEO observed H) some time 

ago that T(1S) -> J/ip + X - 1.1 x 10" 3 . The rate is still 
not sufficiently accounted for |fTo|| , but the most striking 
feature from data is the relatively soft Pj/^ spectrum that 
seemingly peaks below 2 GeV, as shown in Fig. 1(b), 
compared to p max — 4.2 GeV. Perturbative production 
predicts a hard J/tp spectrum [0 (solid curve). Thus, 
this old CLEO result, if confirmed, would be a major 
puzzle for the color octet mechanism. 

Figs. 1(a) and 1(b) together suggest some additional 
mechanism that is responsible for soft J/ip production 
from heavy meson decays. We propose one possibility, 
namely intrinsic charm. IC of hadrons in principle should 
exist Q], and has been suggested to account for charm 
production in deep inelastic scattering Hand J/t/j — > pn 
decay jll]]. It has also been suggested fll^ ] for D and B 
mesons though never pursued. 

For the proton, \p) = \uud) + ^ p uScud \uccud) +■■■, 
the \uccud) component [Lj] is generated by virtual gg — > 
cc interactions (so multi-connected to valence quarks), 
and should scale as a 2 (m 2 )/m 2 relative to the \uud) 
component. Since this higher Fock component arises as 
a quantum fluctuation, ^ p u5cud oc l/(m 2 — M 2 ) where 
M 2 is the invariant mass of the fluctuation. We show 
in Fig. 3(a) the distributions in the \uccud) component. 
One sees that IC carries large momentum fraction Q, in 
contrast to the small x tendency of usual "sea" (extrin- 
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FIG. 3. Intrinsic charm in (a) proton and (b) B meson (b 
and q distributions also shown), with arbitrary normalization. 
Dashed lines in (b) are for the \bccb) component of T(IS'). 

sic) quarks from gluon splitting. Using this feature, there 
is some evidence from data that IC in proton, l^g^J 2 , 
could be at ~ 0.86% level ||. Such an analysis, of course, 
should be done in a more consistent framework E3] and 
incorporate more data to be conclusive. 

For B mesons, one has the differential probability 

dP* k ^(mMl-Eti^) (2) 

dxi ■ ■ ■ dx4 (wo — fh?/xb — rh 2 c jx c — fh 2 /x 5 ) 2 ' 

in the \bccq) Fock component of \B) — ^bq \bq) + 
^becq \bccq) + • • •, since ms, mj, cannot be ignored. In 
fact, the heaviness of guarantees that it still car- 
ries the largest momentum fraction, as can be seen 
from Fig. 3(b). We find (x c ) ~ 0.22 in B, lower than 
(x c ) ~ 0.28 in proton (Fig. 3(a)). 

One has no deep inelastic scattering data off B mesons 
to extract "Fbecgl 2 - However, IvP&gcgl 2 may be no less than 
\^uccud\ 2 , because Jl5| of a larger reduced mass: the B 
meson is more compact than usual hadrons, hence the IC 
amplitude could be dynamically enhanced. Thus, IC in 
B could also be ~ 1%. A heavy quark mass expansion 
study gives IC inp ~ 10~ 3 , which is not inconsistent. 
As stressed in pH}] , this study shows that the cc pair 
should be in the color octet configuration, hence the IC 
component arises from the nonabelian nature of QCD. 

To account for the low pj/^p bump of Fig. 1(a) from the 
\bccq) Fock component, note that B — > J/ipD^ decay 
via Fig. 2(c) is still suppressed by /s, just like Fig. 2(a). 
The spectator decay of Fig. 2(d) is more promising. It 
gives B J/ipDit^ if one assumes factorization, hence 
fits our interpretation of the low pj/^ bump in Fig. 1(a) 
rather well. But can we account for the rate? 

The quark level 3-body spectator decay is very sensi- 
tive to Xb as it scales with the available energy AE to 
the fifth power. Since ~ 0.41 from Fig. 3(b), taking 
Xb ^ 0.41 to 0.6 gives a rate in the 10~ 3 to 5 x 10~ 2 
range, shooting up to even 10% for x b — 0.65. We have 
used AE ~ m b — m c ~ 3.4 GeV for standard b decay 
via \bq) component, and AE ~ X{,to_b — 1-3 GeV via 
\bccq) Fock component of Fig. 2(d). Varying m c ~ 1.3 
to 1.4 GeV does not change this range by much. With 
x c {x c ) ~ 0.22 and xq at its peak near zero, b — > due 



decay leads to the configuration (du)cc(cq), where we 
have indicated color singlet pairings. The remaining cc 
is also color singlet and plausibly evolves into J/tp or r\ c . 
With all quarks having low momenta, there is ample time 
for them to redistribute energy and momentum, includ- 
ing J/ipD + ir~ — > J/ipD Q ir°. It is sensible then that one 
can take the above ~ 10~ 2 factor times the IC fraction 
l^hcojl 2 as a rough estimate. Note that, from duality, 
the "quark level" decay rate can only feed into J/ipD* 
or J/ipDir final state. In this way, we find that a rate 
at few xl0~ 4 level is possible, iflC fraction is not much 
less than 1%. One should check experimentally whether 
our suggested signals are there. 

It is useful to identify other effects where IC of B can 
make important impact. Since the be forms a color singlet 
in \bccq) Fock component, be — > sc, du and Iv annihila- 
tion can proceed without helicity suppression in the vec- 
tor channel, as illustrated in Fig. 2(e). They in general 
complement standard channels hence would not be easy 
to distinguish. However, the energy fluctuation argument 
suggest that the vector be system would have near maxi- 
mal mass, and the accompanying cq would end up in the 
lowest energy-momentum state available, that is, a slow 
moving D*. Indeed, we find that (x c ) + (xq) is right at 
the D* mass. One therefore expects the most significant 
impact to be on B — > D*£v near the zero recoil region, 
where a distortion in spectrum would affect the program 
of \ V c b\ extraction. Thus, if IC in B is truly sizable, there 
would be new systematics to \V c b\ determination, hence 
the importance is beyond QCD. 

Turning to Y(15) — » J/ip + X decay, it is clear that 
IC of T(15) naturally gives soft J/tjj spectrum as given 
in Fig. 1(b). Since the bb in the \bccb) Fock component 
should be color octet ]l6] |, it decays via bb — > g* — > qq, 
gg, as shown in Fig. 2(f) (bb — > gg has an extra i-channel 
contribution), while the accompanying cc is color octet 
and carries minimal energy. Thus, one expects the under- 
lying process T(LS) -> (cc)l%+ 2 "jets", where (cc)^ et 
evolves into slow moving DDX or J/ip(n c )X, in contrast 
with production by leading particle effects from c- 
jets, or from g* — > cc splitting. 

The rate is easier to estimate than the B meson weak 
decay case. We estimate T(IS') — > qccq rate via IC com- 
ponent as, 

T q£c q = x r ee ~ 6 (^) V tecE | 2 x 1.3 keV, (3) 

1 ee V Ci / 

for a single quark flavor q, where 6 is from ratio of color 
and electric charges. Counting 3 (u, d and s) flavors 
and roughly Nc — 3 colors for bb — » gg, one gains an 
additional factor of 6. Understandably, the process of 
Fig. 2(f) is rather fast! We should certainly demand that 
Fqccq < 50% x r^ns), which implies that the IC fraction 
l^bccbl 2 ~ 10 ~ 3 - This is not m conflict with |*h Ec q| 2 ~ 
1% since the T(15) is basically a nonrelativisitic bound 
state hence lacks high frequency components. But it does 
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mean that, if IC is relevant at all, T (IS) -> (cc)^ t + gg 
could easily be 10% (or 5 keV) ofT{\S) rate. 

It is remarkable that the last statement is consis- 
tent with all known facts. Very few hadronic decays of 
T(15) have been reconstructed so far, while the strong 
a| dependence of T(15) — > ggg rate certainly allows 
T(15) -> (cc)^ t + qq ~ 10%. This is also consistent 
with the observed T(15) -> J/V> + X - 1.1 x 1(T 3 , 
since J/^i formation should be just a small fraction of 
(cc)o° t ™ t , while the softness of observed J/ip spectrum is 
also explained. The main part of (cc)g C ™ t evolves into 
£)(*)£)(*) +x where the mesons are slow, accompa- 
nied by two "jets" from the qq and gg with 5 — 6 GeV 
total energy. We note that the ARGUS Collaboration 
has set a bound of @ T(15) -> D*~X < 1.9% for 
Pd* > 0.86 GeV. This is not yet constraining after ad- 
justing for D* fraction, and in particular, one must loosen 
the cut on pn- to be sensitive to IC induced decays. It 
is rather exciting that CLEO would be running on T 
resonances soon jis[| , where we may learn about relative 
abundance of T(15) -» (cc)^ t + 2 "jet" events. Note 
that ~ 1/6 — 1/10 of these "jets" would be charm jets. 

Some further remarks are in order. First, the lump 
below pj/^p < 0.9 GeV cannot be due to feed-down from 
additional cc states, since such spectrum would in general 
extend beyond 1 GeV, as can be seen from data for 
B — > ip' + X and r] c + X. There is no evidence for excess 
above 0.9 GeV. Second, the low pj/^ excess could be due 
to B — > J/ip AN || , but then a large baryon production 
probability is needed near thrchold, which is not more 
plausible. Third, for intrinsic strangeness, one can have 
analgous signatures such as B — > (pD^*> , by making sim- 
ple changes to Figs. 2(a)-(d). While these signals should 
be searched for, they are less distinct since strangeness 
production is less suppressed, and multi-particle final 
states may be preferred since available energy is more 
than 2 GeV. An earlier suggestion ||| of D° -> 4>K Q 
as evidence for intrinsic strangeness in D mesons suffers 
from large rescattering at mo scale. We mention, how- 
ever, the possibility of B — * J/ip<f>K + X, since CLEO 
recently observed |1{| B — > J/ipcpK ~ 10~ 4 . Whether 
it is from intrinsic strangeness or final state rescatter- 
ing, this possibility should be studied further. Fourth, 
analagous signatures for B c mesons and 6-baryons are 
B c -> J/<ipJ/ip(ir), J/ipr) c n and A 6 -> J/V>A c (tt). 

We emphasize that the B — > J/ipD + n~ , J/ipD°ir and 
J/ipD*° signals can be searched for right away at B Fac- 
tories and at the Tevatron with relative easy. The two en- 
vironments complement each other, with excellent ir° de- 
tection at B Factories but larger cross section and higher 
boost at Tevatron, where a large number of J/ip events 
have been recorded. The B — > Dg modes should 

be searched for as control on rescattering diagrams. 

In conclusion, we propose the following study plan. B 
Factory and Tevatron experiments should scan low mo- 



mentum J/ip events for B — > J/ip Dir, J/tp(f>K + X, and 
J/ipAN. If B -> J/ipDir (and perhaps J/^D*) is estab- 
lished above 10 -4 level, while Dg modes are found 
to be far less, one then has a good case for intrinsic charm 
in B. CLEO should run on T(IS') and collect at least 
1 fb _1 data, reconfirm T(15) — * J/ip(r) c ) + X, and search 
for T(15) -> £> ( * } + X. Slow moving J/ip or D^D^ 
in an otherwise jetty event from T(1S*) decay would in- 
dicate, perhaps more unequivocally, the existence of in- 
trinsic charm. The discovery of intrinsic charm in B and 
T mesons would not only add another twist to hadron 
structure, it would have implications on V c b extraction 
from B — > D*lv decay as well. 
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